A phenomenological theory giving the full statistics of the position of fluctuating 
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We propose a phenomenological description for the effect of a weak noise on the position of 
a front described by the Fisher-Kolmogorov-Petrovsky-Piscounov equation or any other travelling 
wave equation in the same class. Our scenario is based on four hypotheses on the relevant mechanism 
for the diffusion of the front. Our parameter-free analytical predictions for the velocity of the front, 
its diffusion constant and higher cumulants of its position agree with numerical simulations. 
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The Fisher Kolmogorov-Petrovsky-Piscounov (FKPP) 
equation 



d t h = d^h + h-h 2 , 



(1) 



describes how a stable phase (h(x, t) = 1 for x — ► — oo) 
invades an unstable phase (h(x, t) = for x — > +oo) 
and how the front between these two phases builds up 
and travels Q. This equation was first introduced in 
a problem of genetics, but equations similar to Q ap- 
pear in much broader contexts like reaction-diffusion 
problems^, optimization!]! , disordered systems^, 
and even particle physics pi I9lll0||. A remarkable example 
is the problem of the high energy scattering of a projectile 
consisting of a small color dipole on a target in the frame- 
work of quantum chromodynamics (QCD): in Ref. Q it 
was recognized that the Balitsky-Kovchegov (BK) equa- 
tion [2j , a mean field equation for high energy scattering 
in QCD, is in the same class as the FKPP equation with 
h being the scattering amplitude, t the rapidity of the 
scattering and x the logarithm of the inverse projectile 
size. 

It is well known|2, El that equations like JIJ have a 
family of travelling wave solution of the form h{x, t) — 
h(z) with z = x — vt. There is a relation between the 
exponential decay of each solution (h(z) ~ exp(— jz) for 
large z) and its velocity: v = v(j). For example, 17(7) = 
7+I/7 for the FKPP equation QJ. Other front equations 
would give different expressions of 11(7) . See, for example, 
section |IV| or Refs. [ll [3 . 

If one starts with a steep enough initial condition, the 
front converges to the travelling wave with the minimal 
velocity. Therefore 

^deterministic = minwf/y) = 11(70) where 1/(70) = 0, 
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^deterministic (^) ~ AzG . 



(2) 



(The multiplicative factor z in ^deterministic is present only 
for this slowest moving solution.) 

There is a large class (the FKPP class) of equations de- 
scribing the propagation of a front into an unstable state 



which select the minimal velocity, as described by 
(There exist also equations, called "pushed" or "type II" 
for which the velocity selected by the front is not the 
slowest one. The properties of these fronts are quite 
different 0, 0] from the properties of QJ, and we will 
not consider them in the present paper.) 

Deterministic front equations such as QJ usually oc- 
cur as the limit of a stochastic reaction-diffusion model 
[l5j when the number of particles (or bacterias, or reac- 
tants) involved becomes infinite. In a physical situation, 
all numbers remain finite and a small noise term should 
be added to (JTJ to represent the fluctuations at the mi- 
croscopic scale. One might write, for instance|l6j. 



d t h = dlh + h-h 2 + y/h(l - h) /N r?(x, t), (3) 

where rj(x,t) is a normalized Gaussian white noise and 
N is the number of particles involved. 

The effect of such a noise is to make the shape of the 
travelling wave fluctuate in timeQ. It affects also its 
velocity and makes the front diffuse 0,0,0- 

For a chemical problem, N might be of the order of 
the Avogadro number and one could think that such a 
small noise term should give small corrections, of order 
l/ViV, to the shape and position of the front. However, 
because the front motion is extremely sensitive to small 
fluctuations in the region where h ~ 1/N, this is not the 
case. In presence of noise as in (jHJ, the front has an ex- 
ponential decay if h(x,t) 1/N, but it vanishes much 
faster than this exponential in the region where h(x, t) is 
of order 1/N Q. (This is obvious in a particle model, as 
there cannot be less than one particle at a given place.) 
As an approximation to understand the effect of the mi- 
croscopic stochastic details of the system, it has been 
suggested to replace the noise term by a deterministic 
cutoff which makes the front vanish very quickly when 
h ~ 1/N For instance, for the FKPP equation JTJ, 
one way of introducing the cutoff is 



d t h = d 2 x h + (h - h 2 )a(Nh), 
with a(r) = 1 for r > 1 and a(r) — > for r 



(4) 



In the presence of such a cutoff, the velocity and the 
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shape 0) become, for any equation in the FKPP class, 

hcntos(z) ~ A— sin f ^) e~ 7 ° 2 where L = —log TV, 

(5a) 

(5b) 



^cutoff ~v(j ) 
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(The shape (|5a() is valid only in the linear region, where 
h is small enough for the nonlinear term h 2 to be neg- 
ligible but still larger than 1/N. Note that for z <C L, 
the shape coincides with 0. A way to interpret the sine 
is to say that the front moves slower than the minimal 
velocity ^deterministic = v(lo) and that the decay rate be- 
comes complex: 7 = 70 ± iir/L. Then, the expression of 
^cutoff results from an expansion of 1^(7) for large L.) 

The prediction does not depend on the details of 
the microscopic model. It only depends on the deter- 
ministic equation and on the existence of a microscopic 
scale. This cutoff picture is also present in the mean field 
QCD context in [18j, where it was introduced to avoid 
unitarity violating effects in the BK equation at interme- 
diate stages of rapidity evolution. In this context, N is 
i/ a QCD where oqcd is the strong coupling constant. 

Extensive numerical simulations of noisy fronts have 
been performed over the years [3. Il7j. and the large cor- 
rection l|5bjl to the velocity found in the cutoff picture 
seems to give the correct leading correction to the veloc- 
ity of noisy fronts. (See [l^ for rigorous bounds.) Being 
a deterministic approximation, the cutoff theory gives 
however no prediction for the diffusion constant of the 
front. 

In the present paper, we develop a phenomenological 
description which leads to a prediction for this diffusion 
constant. This description tries to capture the rare rele- 
vant events which give the dominant contribution to the 
fluctuations in the position of the front. The prediction 
is that the full statistics of the front position in the noisy 
model depends only on the amplitude 1/N of the noise 
at the microscopic scale and on 11(7), a property of the 
deterministic equation. For large N, all the other details 
of the underlying microscopic model do not contribute to 
the leading order. Our description leads to the following 
prediction for the velocity and for the diffusion constant 
of the front for large N: 



V - Ucutoff = 7T 7 U (70) 



D = ^ 7o V'(7o) 



3 log log N 
70 log 3 N 

7T 2 /3 

7o log N 



(6a) 



Actually, our phenomenological approach also gives a 
prediction to the leading order for all the cumulants of 
the position of the front. For n > 2, 



-th cumulant] 



^7 V'(7o) 



IS lc >g 3 N 



(6b) 



where £(n) = Ylk>i ^ 

The 1 / log 3 N dependence of the diffusion constant was 
already observed in numerical simulations |17|. In the 
QCD context, it was proposed in [23 to identify the full 
QCD problem with a stochastic evolution, such as 10, 
and the dependence of the diffusion constant was used 
to suggest a new scaling law for QCD hard scattering at, 
perhaps, ultrahigh energies. 

We do not have, at present, a mathematical proof of 
the results (|6af> and (|6b|l . Rather, we believe that we have 
identified the main effects contributing to the diffusion of 
the front. We present our scenario in Sec. ITll where we 
state a set of four hypotheses from which the results © 
follow. We give arguments to support these hypotheses 
in sections |III Al to IIII Dl Finally, to check our claims, 
we present numerical simulations in section IIVI for the 
five first cumulants of the position of the front. These 
simulations match very well the predictions 0. 



II. THE PICTURE AND ITS QUANTITATIVE 
CONSEQUENCES 

To simplify the discussion, we consider, in this section, 
more specifically a microscopic particle model rather 
than a continuous stochastic model such as 0. This 
is merely a convenience to make our point clearer, but 
the discussion below could be rephrased for other models 
in the stochastic FKPP class. 

We consider models where particles diffuse on the line 
and, occasionally, duplicate. If one considers, for h(x,t), 
the density of particles or, alternatively, the number of 
particles on the right of a;, it is clear that it is not yet 
described by a front equation, because it grows exponen- 
tially fast with time; one needs to introduce a satura- 
tion rule. For instance: I) keep the number of particles 
fixed by removing the leftmost particles if necessary; or 
2) remove all the particles which are at a distance larger 
than L behind the rightmost particle; or 3) limit the 
density by allowing, with a small probability, that two 
particles meeting recombine into one single particle 



A. A scenario for the propagation of the front 

The main picture of our phenomenological description 
is the following. The evolution of the front is essentially 
deterministic, and its typical shape and velocity are given 
by Eq. Q. But from time to time, a fluctuation sends 
a small number of particles at some distance S ahead of 
the front. At first, the position of the front, determined 
by where most of the particles are, is only modified by a 
negligible amount of order I /N by this fluctuation. How- 
ever, as the system relaxes, the number of wandering par- 
ticles grows exponentially and they start contributing to 
the position of the front. Meanwhile the bulk catches 
up and absorbs the wandering particles and their many 
offsprings; finally, the front relaxes back to its typical 
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shape J5aJ. The net effect of a fluctuation is therefore 
to shift the position of the front by some amount R(5) 
which depends, obviously, on the size S of the fluctua- 
tion. A useful quantity to characterize the fluctuations 
is the width of the front, ft can easily be defined as the 
distance between the leading particle (where h w 1/N) 
and some position in the bulk of the front, for instance, 
where h — 0.5. (Changing this reference point would 
change the width by a finite amount, independent of N). 
This width is typically of order L, where L is given by 
the cutoff theory l|5al) . During a fluctuation that sends 
particles at a distance 5 ahead of the front, the width 
of the front increases quickly to L + 5, and then relaxes 
slowly back to L. 

We emphasize that, in this scenario, the effect of noise 
is so weak that, most of time, it can be ignored and 
the cutoff theory describes accurately the evolution of 
the front. It is only occasionally, when a rare sequence 
of random microscopic events sends some particles well 
ahead of the front that the cutoff theory is no longer 
valid. The way this fluctuation relaxes is, however, well 
described by the deterministic cutoff theory. 

We shall encode this scenario in the following quanti- 
tative assumptions: 

1. If we write the instantaneous fluctuating width of 
the front as L + 6, then the probability distribution 
function for 6 is given by 



p(6)dS = de'^dd, 



(7) 



where G\ is some constant. Note that we assume 
this form only over some relevant range of values: 
6 large enough (compared to I) but much smaller 
than L (typically of order logi). Fluctuations 
where 6 is "too small" are frequent but do not con- 
tribute much to the front position. Fluctuations 
where 5 is "too large" are so rare that we do not 
need to take them into account. Only for "moder- 
ate" values of S do we assume the above exponential 
probability distribution function. 

2. The long term effect of a fluctuation of size 6 (as- 
suming that there are no other fluctuations in- 
between) is a shift of the front position by the quan- 
tity 

P 7o<5' 



R(S) = -log(l + C 2 — 
7o V L 

where C2 is another constant. 



(8) 



3. The fluctuations of the position of the front are 
dominated by large and rare fluctuations of the 
shape of the front. We assume that they are rare 
enough that a given relevant fluctuation has enough 
time to relax before another one occurs. 

From these three hypotheses alone, one can derive our 
results (0 up to a single multiplicative constant. This 
constant can be determined with the help of a fourth 
hypothesis: 



4. For the aim of computing the first correction to the 
front velocity obtained in the cutoff theory (J5J , one 
can simply use the expression <|5b(l with L replaced 
by L e ff where 



f 3 

L cS = — log TV + — log log iV - 

7o 7o 



(9) 



It is important to appreciate that the average or 
typical width of the front is still L and not L c ff- 
The latter quantity is just what should be used in 
l|5b(l to give the correct velocity. 

B. How ||6j) follows from these hypotheses 

We are now going to see how the results © follow from 
these four hypotheses. 

First, we argue that the probability to observe a fluc- 
tuation of size S during a time interval At can be written 
as p(S)d5 At/r, where p(S) is the distribution J7J) of the 
increase of the width of the front and where r is some 
typical time characterizing the rate at which these fluctu- 
ations occur. Indeed, during a fluctuation of a given size, 
the width of the front increases to that size and then re- 
laxes back. For a large 8, observing a front of size L + S is 
very rare, but, when it happens, the most probable is that 
one is observing the maximum expansion of a fluctuation 
with a size close to 5; the contribution from fluctuations 
of sizes significantly larger than S is negligible as they are 
much less likely. 

Second, as a fluctuation builds up at the very tip of 
the front where the saturation rule (see beginning of sec- 
tion [nj can be neglected, we argue that the typical time 
r introduced in the previous paragraph and the time it 
takes to build a fluctuation of a given size do not depend 
on N. (However, the relaxation time of a fluctuation de- 
pends on N as the bulk of the front is involved in the 
relaxation.) 

Let X t be the position of the front, i5o the minimal size 
of a fluctuation giving a relevant contribution to the po- 
sition of the front and At a time much smaller than the 
time between two relevant fluctuations, but much larger 
than the time it takes to build up such a fluctuation and 
have it relax. (This is authorized by the third hypothe- 
sis.) We have 

Xt+At — 

(X t + v cutoS At + R(S) proba. ^p(5)dS for S 



X t + v cutoS At 



proba. 1-M f™p(S)d6 



(Note that =^ p(S)dS is the probability of observing 
a relevant fluctuation during the time At. By definition 
of At, this is much smaller than I.) 

One can then compute the average, denoted (•), of 
exp(XX t +At)- One gets, for A small enough, 



d t log ( 



) = Aw cutofF +^ / p{5) 



dS. (10) 
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Expanding in powers of A, one recognizes on the left- 
hand-side the cumulants of X t . Therefore, one gets 



V - Vcutoff = - / i>( '>)/?('>) M. 

[n-th cumulant] = 1 ; pfflg , ffl # for „ > 2 



(11) 



At this point, one can notice from the expressions of p(S) 
and R(S) that the values of S such that e 705 3> L 3 have a 
negligible contribution to the integrals giving the velocity 
and the cumulants. Thus appears naturally a 5 max = 
(3/7o) log L which is exactly the effective correction to 
the width of the front appearing in |J3J| . 

The integrals in Eq. Ijllfl can be evaluated, and one 
gets 



p(S)R n (S) dd 



C1C2 



log" 1 



1 



dx, 



with x = (L 3 /C 2 ) exp(-7 (5) 
gives log(L 3 /C 2 ). For n > 2 
to 00 (the correction is at most of order 1/L 6 ) and one 
recognizes n!£(n). Finally, 



(12) 

For n = 1, this integral 
one can integrate from 



V - Wcutoff 



C*iC 2 31og£ 
T70 70 L 3 
[n-th cumulant] C\Gi n\Q(n) 



(13) 



t 



76 



L 3 



Everything is determined up to one numerical con- 
stant C1C2/V. As the fourth hypothesis gives the ve- 
locity, one can easily determine that constant and re- 
cover (0. 

All the cumulants (except the first one) are of the same 
order of magnitude, as the fluctuations are due to rare 
big events. 



B. Second hypothesis 

To obtain (JSJ), we need to compute the response of 
the deterministic model with a cutoff jlj to a fluctuation 
at the tip of the front. This is a purely deterministic 
problem: starting with a fluctuation (i.e. a configuration 
slightly different from the stationary shape), we let the 
system evolve with a cutoff and relax back to its station- 
ary shape H5a[) . and we would like to compute the shift 
in position due to this fluctuation. 

Although the evolution is purely deterministic, the 
problem remains a difficult one. For simplicity, we discuss 
here the case of the FKPP equation QJ. The extension 
to other travelling wave equations in the FKPP class is 
straightforward. 

There are two non-linearities in (£!}: one is the ~h 2 
term, which is important when h is of order 1, and the 
other one is the cutoff term a(Nh), which is important 
when h is of order 1/N. Between these two points, there 
is a large length of order L = log N where one can neglect 
both non-linearities. This means that, for all practical 
purpose, one can simply use the linearized version of the 
FKPP equation for the whole front except for two small 
regions with a size of order 1 at both ends of the front. 

Let X t be the position of the front, and L t its length. 
There are many equivalent ways of defining precisely 
these quantities; for instance we can take X t such that 
h(X t ,t) = 10~ 5 and L t such that h(X t + L t ,t) = ±. We 
expect that X t — w cu toff t and L t — L, which are quantities 
of order 1, have a relaxation time of order L 2 , as for the 
shape of the front. 

For Xt < x < X t + L t , the problem is linear: 



d t h = dlh + h. 



Using the Ansatz 
h(x,t) = L t G 



Li 



Xt t_ 

' Li 1 



_ ) p -0-fcutofft) 



(14) 



(15) 



III. ARGUMENTS TO SUPPORT THE 
HYPOTHESES 

A. First hypothesis 

This first hypothesis is not very surprising if one con- 
siders that exp(— 7o<5) is the natural decay rate of the 
deterministic equation. A more quantitative way to un- 
derstand Q is that building up a fluctuation is an effect 
which is very localized at the tip of the front, where satu- 
ration effects can be neglected. We present in appendix lAl 
a calculation using this property. 

Moreover, numerical simulations |2l| of that proba- 
bility distribution function give evidence that for large 
enough TV, the decay is exponential with the rate 70 as 
in 0. 



with ^cutoff = 2 - (see IgbJ for ^(7) = 7 + I/7) 
and keeping only the dominant terms in L, the func- 
tion G(u,t) evolves according to 

^ 2r ' < ir 2 G, 



d 2 G 



d T G 

with the boundary conditions 



G(0,t)«0, G(1,t)«0. 



(16) 



(17) 



(More precisely, G(0,t) and G(1,t) would be non zero 
only at the next order in a 1/L expansion.) 

The problem reduces to a diffusion problem with ab- 
sorbing boundary conditions. The stationary configura- 
tion is the sine shape l|5all . as expected. 

If, at time t = the shape is different from this sta- 
tionary configuration, it will relax back to it in the long 
time limit up to a multiplicative constant: 



j2 

G(y, 00) = — sin(7ry). 

7T 



(18) 
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As the stationary shape for h(x, t) must be of the form 
given by l|5a[) . we obtain, using (|15fl that the final shift 
in position is given by 



R(S) = lim (X t - Vcutost) 



(19) 



To compute the value of B, one simply needs to project 
the initial condition on the sine shape: 

B = Ae R{&) = 2vr f dy Bm(wy)G(y,0). (20) 
Jo 

We now proceed to use this expression for the pertur- 
bations we are interested in: perturbations localized near 
the cutoff. 

We do not have a full information on the initial con- 
dition h(x,0) or, equivalently, G(y,0). However, as we 
expect a perturbation to grow at the very tip of the 
front, we expect that h(x,Q) is identical to its station- 
ary shape, except in a region of size of order Aa; ~ 1 on 
its tip. On the scale we consider, this means that G(y, 0) 
is perturbed over a region of size Ay w 1/L. In other 
words: 



G(y,0) = A 



1 



sin(Try) + p(l - y) 



(21) 



where the perturbation p{y') is non-zero only for y' 
1 — y of order 1/L. Therefore, from J3UJ, 



,R(S) 



1 + 2tt 



dy' ny'p(y'), 



(22) 



where b is a number of order 1 representing the extent 
over which a perturbation initially affects the shape of 
the front, (p(y') w if y' > b/L.) 

The precise shape of p(y') is not known, but its am- 
plitude can be easily understood in a stochastic particle 
model: if some particles are sent at a distance 8 <C L 
ahead of the front, h(x,t) increases by 1/N at position 
x = X t + L + 5. Because of the exponential factor in (|15|l . 
this translates to an increase of order p(y') « exp(6)/L 
for the reduced shape G(y,r). Combining everything, 
one finally gets 



Ms) 



l + C 2 



L 3 ' 



(23) 



where C2 is some number of order 1 which depends on the 
precise shape p(y'). Expression (|23|) is just our second 
hypothesis, up to factors 70 which can be put back by 
dimensional analysis. 

One consequence of the argument above is that Ci is of 
order 1 compared to L. However, it gives no information 
about the dependence of C2 on 8 or on the shape of the 
fluctuation. We think that if C2 depends on 8, it is a 
weak dependence that we can ignore. A simple situation 
where this can be checked is when 8 is large: if a particle 
jumps sufficiently far ahead, it will start a front of its own 



that will completely replace the original front. For such a 
front, it is well known [121I22T ] that the position for large t 
is given at first (while the cutoff is not relevant) by 8 + 
It — I logf. When the velocity 2— ^ matches i> C utoff, that 
is at a time to « L 2 , a crossover occurs and the position 
becomes R(S)+v cuto s t. Matching the two expressions for 
the position at time t = to, one obtains R(8) ~ 8 — logi 3 , 
as predicted by (JHJ. This indicates that, at least for large 
8, the number C2 has no 8 dependence. 



C. Third hypothesis 

From section |n] and Eq. JHJ) , the size 8 of the fluctua- 
tions that contribute significantly to the diffusion of the 
front are such that exp(7o<5) ~ L 3 . From J7J, the typical 
time between two such fluctuations is therefore L 3 . On 
the other hand, from section IIII Bl the relaxation time 
of a fluctuation is of order L 2 . It is therefore safe to as- 
sume that a relevant fluctuation has enough time to relax 
before another one occurs. 



D. Fourth hypothesis 

The fourth hypothesis states that, to compute the shift 
in velocity, one should use a front width L e g that is larger 
than what is predicted by the cutoff theory by an amount 
^■loglogiV. The hypothesis is plausible as this length is 
precisely the distance 8 at which the relevant fluctuations 
occur: the main effect of the fluctuations would then be 
to increase the effective width of the front that enters the 
cutoff theory (J5J. We present in appendix iBl a simplified 
model to support this claim. 

Remarkably, the front width L e ff emerges naturally in 
the QCD context . 



IV. NUMERICAL SIMULATIONS 

We consider here a reaction-diffusion model with sat- 
uration which was introduced in |_13j as a toy model for 
high energy scattering in QCD. Particles are evolving 
in discrete time on a one-dimensional lattice. At each 
timestep, a particle may jump to the nearest position on 
the left or on the right with respective probabilities pi and 
p r , and may divide into two particles with probability A. 
We also impose that each of the n(x, t) particles piled up 
at x at time t may die with probability Xn(x,t)/N. 

Between times t and t + 1, ni(x,t) particles out of 
n(x,t) move to the left and n r (x,t) move to the right. 
Furthermore, n+(x, t) particles are replaced by their two 
offsprings at x, and n_(a;, t) particles disappear. Hence 
the total variation in the number of particles on site x 
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model I 
model II 

asymptotic prediction (9) 
prediction with sub-leading terms (31) 



N 



FIG. 1: Measured L e s, defined by 1301 . from which we have subtracted the width L in the cutoff theory, as a function of N. 
The dashed line represents the leading terms 3 log log iV/70 , see @ . The subleading terms 13111 of the plain line have been 
determined by a fit. 



reads 



For the purpose of our numerical study, we set 



n(x, t + 1) — n(x, t) = —rii(x, t) — n r (x, t) — n_(x, t) 

+ n+(x,t) + ni(x + l,t) + n r (x - l,t). (24a) 

The numbers describing a timestep at position x have a 
multinomial distribution: 



P({ni,n r ,n + ,ri-}) = 



1! 



rPj Pr 



n;!n r !n+!n_!An! 
A"+ (An/iV) n - (1 -pi -p r -X-\n/N) An , (24b) 

where An = n — ni — n r — n + — n_ , and all quantities in 
the previous equation are understood at site x and time 
t. The mean evolution of u = n/N in one step of time 
reads 

{u(x, t+l)\{u(x, t)}) = u(x, t) +pi[u(x + l, t) — u(x, t)] 
+p r [u(x-l, t)-u(x, t)]+Xu(x, t)[l-u(x, t)}. (25) 

When N is infinitely large, one can replace the u's in l|25|) 
by their averages. One obtains then a deterministic front 
equation in the FKPP class with 

«( 7 ) = - log (1 + A + j3/(e~ 7 - 1) + Pr (e-> - 1)) , (26) 
7 

and 70 is defined by 1/(70) = 0, see @- 



p /= p r = 0.1 and A = 0.2. 

From H26f) . this choice leads to 

70 = 1.3521 ■ ■ • , u(7o) = 0.25538 • 
u"( 7o ) = 0.16773- •• . 



(27) 



(28) 



Predictions for all cumulants of the position of the front 
are obtained by replacing the values of these parameters 
in ©. 

Technically, in order to be able to go to very large 
values of N, we replace the full stochastic model by its 
deterministic mean field approximation u — > (it), where 
(u) is given by Eq. 1251) . in all bins in which the num- 
ber of particles is larger than 10 3 (that is, in the bulk of 
the front). Whenever the number of particles is smaller, 
we use the full stochastic evolution (J22J. We add an ap- 
propriate boundary condition on the interface between 
the bins described by the deterministic equation and the 
bins described by the stochastic equation so that the 
flux of particles is conserved |23j. This setup will be 
called "model I" . Eventually, we shall use the mean field 
approximation everywhere except in the rightmost bin 
(model II): at each time step, a new bin is filled imme- 
diately on the right of the rightmost nonempty site with 
a number of particles given by a Poisson law of average 
9 = N(u(x,t+l)\{u(x,t)}). In the context of a slightly 
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FIG. 2: From top to bottom, the correction to the velocity given by the cutoff theory and the cumulants of orders 2 to 5 of the 
position of the front in the stochastic model. The numerical data are compared to our parameter-free analytical predictions JBJ , 
represented by the dashed line. The subleading terms of the plain lines are numerically the same as in figure no further fit 
has been performed for the present figure. 



different model in the same universality class this 
last approximation was shown numerically to give in- 
distinguishable results from those obtained with the full 
stochastic version of the model, as far as the front veloc- 
ity and its diffusion constant were concerned. We shall 
confirm this observation here. 

We define the position of the front at time t by 
00 

X t = Y,u(x,t). (29) 

.T = 

We start at time t = from the initial condition u(x, 0) = 
1 for x < and u(x,0) = for x > 0. We evolve it up 
to time t = log 2 N to get rid of subasymptotic effects 
related to the building of the asymptotic shape of the 
front, and we measure the mean velocity between times 
log 2 N and 16 x log 2 N. For model I (many stochastic 
bins), we average the results over 10 4 such realizations. 
For model II (only one stochastic bin), we generate 10 5 
such realizations for N < 10 50 and 10 4 realizations for 
N > 10 50 . In all our simulations, models I and II give 
numerically indistinguishable results for the values of N 
where both models were simulated, as can be seen on the 
figures (results for model I are represented by a circle and 
for model II by a cross). 

First, we check that the effective width of the front 
is L e g given by Eq. We extract the latter from the 



measured mean velocity v using the formula 

^xlwr^- (30) 

V 2«7o) - v) 

We subtract from L e g the width of the front obtained 
in the cutoff theory L = (\ogN)/j n , and compare the 
numerical result with the analytical formula 

31og(logAQ >g(logiV) ( . 

L cS ~L = hc + d — — — . (31) 

7o log N 

The first term in the r.h.s. is suggested by our fourth 
assumption (see Eq. ©). We have added two sublead- 
ing terms which go beyond our theory: a constant term, 
and a term that vanishes at large N. The latter are nat- 
urally expected to be the next terms in the asymptotic 
expansion for large N. We include them in this numer- 
ical analysis because in the range of N in which we are 
able to perform our numerical simulations, they may still 
bring a significant contribution. 

We fit l|31l) to the numerical data obtained in the 
framework of model II, restricting ourselves to values of 
N larger than 10 30 . In the fit, each data point is weighted 
by the statistical dispersion of its value in our sample of 
data. We obtain a determination of the values of the free 
parameters c = —4.26 ± 0.01 and d — 5.12 ± 0.27, with a 
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FIG. 3: The ratio of cumulants 2 to 5 divided by the correction to the velocity to the power 3/2. The dashed lines are the 
analytical prediction assuming only the cutoff theory (15 bl for the velocity and the prediction JSJ for the cumulants. 



good quality of the fit (x 2 /d.o.f ~ 1.15). The numerical 
data together with the theoretical predictions are shown 
in figure ^ We see a clear convergence of the data to the 
predicted asymptotics at large N (dotted line in the fig- 
ure), but subleading corrections that we have accounted 
for phenomenologically here are sizable over the whole 
range of N. 

We now turn to the higher order cumulants. Our nu- 
merical data is shown in figure |3 together with the ana- 
lytical predictions obtained from JfJJl (dotted lines in the 
figure). We see that the numerical simulations get very 
close to the analytical predictions at large N. However, 
like in the case of L c g, higher order corrections are pre- 
sumably still important for the lowest values of N dis- 
played on the plot. 

We try to account for these corrections by replacing 
the factor (logiV)/7o = L in the denominator of the ex- 
pression for the cumulants in Eq. JBJ by the Ansatz for 
L c ff given in l|31|l . without retuning the parameters. The 
results are shown in figure [21 (full lines), and are in excel- 
lent agreement with the numerical data over the whole 
range of N. We could also have refitted the parameters 
c and d for each cumulant separately, as, a priori, they 
are not predicted by our theory. We observe that this is 
not required by our data. 

This last observation suggests that all the cumulants 
can be computed, with a good accuracy, with the effec- 
tive width L e ff as the only parameter. We check this on 
figure |3 which represents the ratio of the n-th cumu- 
lant (divided by time) by the correction to the velocity 



^deterministic — V to the power 3/2. If one supposes that 
the correction to the velocity varies like l/L^ s and the 
cumulants like 1/L^ S for some effective width L c g, this 
width disappears from the ratio plotted and one can com- 
pare the numerical results to our analytical prediction 
with no free parameter or unknown subleading terms. 
Within statistical error, the data seem to agree for N 
large enough with our prediction, suggesting that, in- 
deed, all the cumulants can be described with a good 
accuracy with only the effective width Left- 
Simulations, not shown here, for the model introduced 
in support also our predictions ©. 



V. CONCLUSION 

The main idea that we have put forward in the present 
work is that all the fluctuations of the front position, 
and in particular the diffusion constant, are dominated 
by large but rare fluctuations at the tip of the front. 

Under some more precise assumptions (hypotheses of 
section [HJ on these fluctuations, we were able to obtain 
explicit expressions JBJ of the cumulants of the position 
of the front. We checked these predictions in our numer- 
ical simulations of section IIVI In section IIIII we gave 
some arguments in support of the four hypotheses of sec- 
tion [H] None of these arguments can be regarded as a 
mathematical derivation, and we can imagine that some 
details, such as the precise shape of the distribution of 
fluctuations (7J or the explicit expression (fHJ , could be 
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slightly modified by a more precise analysis. We believe 
however, given the good agreement of the predictions © 
with the numerical simulations, that our picture is very 
close, if not identical, to the actual behavior of the front 
for large values of N. 

To conclude, we would like to point out the remarkable 
similarity between the predictions 10 and the exact re- 
sults obtained recently [2^ in the context of directed poly- 
mers. Basically, the results of [24| are the same, mutatis 
mutandis as our present results ©, for all the cumulants. 
The only significant change is that the 3 log log N for the 
velocity and the 1 / log 3 N dependence for all the cumu- 
lants in © corresponds, in 24], to a log log N for the 
velocity and a 1/logiV for all the cumulants [24. Eq. (23) 
with L = logiV. The term L + logL in the velocity corre- 
sponds to ^cutoff, see Eq. (28)]. What is interesting is that 
our scenario of section^for FKPP fronts applies also for 
the system studied in [2J]: Indeed, the fluctuations of the 
position are mainly due to the rare big events taking place 
at the tip of the "front" , [24l last paragraph before con- 
clusion] , the position of the rightmost particle is given by 

H E q- ( 32 ) with 6 = -logg and X t = log£ t ], the 
effect of a large fluctuation can be written as JSJ) with the 
L 3 term replaced by L [24l the log of (34) can be written 
as X t +i —X t = L + logL + R(S)], relevant fluctuation (of 
size log L instead of 3 log L) appears every L timesteps 
(instead of every I? timesteps) and the relaxation time is 

1 instead of L 2 . This similarity may add a further piece 
of evidence for our results. 

This work was partially supported by the US Depart- 
ment of Energy. 



where rj is a Gaussian number of variance 2, and dupli- 
cates with a probability dt. If it duplicates, the probabil- 
ity P t +dt{x) is the probability that the offsprings of both 
particles are on the left of x. As the particles have un- 
corrected motion, this is the product of the probabilities 
for each offspring. Finally, one eets[25| 



Pt+dt{x) = {P t {x - vVdt)(l - dt) + P 2 {x - r]Vdt)dt) , 

where the average is on rj. After simplification, 

d t P = d 2 x P-P + P 2 . (Al) 

One notices that 1 — P t (x) is solution of the deterministic 
FKPP equation Therefore, for large t and x, 0.0, 



l-Pt(x) 



ze z 4 * for z = x — 2t + - \ogt. 



Let Qt(x) be the probability that there are no particles 
on the right of x when the initial condition is a given 
density of particles po(x). Using the fact that all the 
particles are independent, one gets easily 



Qt(x) 



exp 



dy Po(y)(l - Pt(x - y) 



(A2) 



Po(y) needs to reproduce the shape of the front seen 
from the tip. Starting from l|5a(l . we write p(y) — 
Nhc U toft(L + y) and take the large N limit. One gets 
Po(y) = ~yexp(-y) for y < and p (y) = for y > 0. 
Evaluating the integral in (|A2|) . one gets, for large t and 
x - 2t < y/i, 



APPENDIX A: LIMIT N ^ oo 

In this appendix, we try to provide an argument for the 
exponential decay J7J) of the distribution for the width of 
the front. To this aim, we consider a very simple model of 
reaction-diffusion: particles diffuse on the line and dur- 
ing each time interval dt, each particle duplicates with 
a probability dt. The motions of all the particles are 
uncorrelated. 

If one added a saturation rule as described at the be- 
ginning of sectional the density of particles (or the num- 
ber of particles on the right of x, depending on the pre- 
cise saturation rule) would be described by a stochastic 
FKPP equation. However, the saturation affects only the 
motion of particles in the bulk of the front, where the 
density is high. As the fluctuations develop in the low 
density region, it is reasonable to assume that the distri- 
bution of the size of the fluctuations are well described 
by the model without any saturation. 

For this model without saturation, let Pt(x) the prob- 
ability that, at time t, no particles are present on the 
right of x given that, at t — 0, there is a single particle at 
the origin: Pq(x) = 9(x). During the first "time step" dt, 
the only particle in the system moves by a quantity r\\fdi 



t(x) « exp 



_ Ce -(x-2t) 



(A3) 



(Notice how the (3/2)logi factor canceled out). The 
probability distribution function of the rightmost particle 
is clearly d x Q t (x). We see that in this stochastic model, 
the front moves at a deterministic velocity equal to 2 and 
that the position of the rightmost particle around the po- 
sition of the front is given by a Gumbel distribution. 

The distribution (|A3|) gives our first hypothesis (Q for 
large fluctuations (5 = x — 2t 3> 1). Our attempts to 
check numerically (|A3I) by simulating fronts with a large 
but finite number of particles confirmed this exponential 
decay for large 6, but showed some discrepancy for 6 < 0, 
which we do not understand. This, however, does not 
affect the hypothesis Q. 



APPENDIX B: MOVING WALL 

We consider again the reaction-diffusion model intro- 
duced in appendix lAl As we said, one needs to add a 
saturation effect to obtain a propagating front equation 
for the density, but doing so introduces correlations in 
the motions of the particles that make the model hard 
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to solve. In this appendix, we introduce an approximate 
way of adding a saturation effect which does not intro- 
duce any such correlation. 

In a real front, the tip is subject to huge fluctuations 
happening on short time scales. On the other hand, the 
bulk of the front moves smoothly and adjusts very slowly 
to the fluctuations happening at the tip. Therefore, we 
believe that, for times not too large, it is a reasonable 
approximation to assume that the bulk of the front moves 
at a constant velocity. 

To implement this idea, our model is the following: a 
wall starting at the origin is moving to the right at a 
constant velocity v. Particles are present on the right of 
the wall. The particles are evolving as in appendix ^ 
except that whenever a particle crosses the wall, it is 
removed. 

We first consider a single particle starting at a dis- 
tance z of the wall. After a time t, either all the off- 
springs of this particle have been caught by the wall, or 
some have survived. We want to compute the proba- 
bility E t (z) that all the particles have been caught at 
time t. The original particle, after a time dt, is at a dis- 
tance z — vdt + rjVdt from the wall, and it might have 
duplicated with probability dt. Using the same method 
as in appendix lAl one gets 



d t E t = dlE t - vd z E t -Et + El 



(Bl) 



with the conditions 



E (z) = for z > and E t {z) = 1 for z < 0. (B2) 
In the long time limit, E t (z) converges to the station- 



ary solution E 00 (z), and one recognizes that h(z) — 
1 — Eoc^z) is the stationary solution of the FKPP equa- 
tion JJJ if z = x — vt. In other words, 1 — £' 00 (— z) is 
the shape of a travelling front. As this shape reaches 
for z = 0, it must be a front with a sine arch and a 
velocity v smaller than 2, as in ©. So if v < 2, the 
probability 1 — i? QO (— z) is the shape of the front with a 
cutoff: 



1 - E(z) ~ Le 



L sin {j^y) eZ wnere v ~ ^ — j-^ ■ 

(B3) 

(The extra factor e L comes from the fact that z — 
is the tip of the front in (|B3() while it is the bulk of the 
front in 15a|) . If v > 2, all the particles eventually die 
and E t (z) converges to 1.) 

If one starts with a density p{z) of particles at time t = 
0, the probability E£ that everybody dies is given, simi- 
larly to HA2jl. by 



Ef = exp 



dz p(z)(l-Et(zj) 



(B4) 



We consider, as an initial condition, the situation in the 
real front with p(z) = Nh(z) ~ NL sin(7rz/i) exp(— z) 
for z < L as in J5aj. One gets, for long times, 



E% -> exp [-CNL 3 e- L ] 
We see that the system survives if 



(B5) 



NL 3 e- L > 1 or L<\ogN + 3log\ogN, (B6) 
which, given 70 = 1, is exactly (JSJ. 
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